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Abstract 

We present a new connection between the Hele-Shaw flow, also known as two- 
dimensional (2D) Laplacian growth, and the theory of holomorphic discs with 
boundary contained in a totally real submanifold. Using this we prove short time 
existence and uniqueness of the Hele-Shaw flow with varying permeability both 
when starting from a smooth Jordan domain and when starting from a point. Ap- 
plying the same ideas we prove that the moduli space of smooth quadrature do- 
mains is a smooth manifold whose dimension we also calculate, and we give a 
local existence theorem for the inverse potential problem in the plane. 



1 Introduction 

The Hele-Shaw flow is a model for describing the propagation of fluid in a Hele- 
Shaw cell. Such a cell consists of two parallel plates separated by a small gap, and 
the fluid is confined to the narrow space between them. This moving boundary 
model has been intensely studied for over a century, and is a paradigm for under- 
standing more complicated systems such as the flow of water in a porous media, 
the melting of ice and models of tumour growth. 
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Figure 1 : Cross section of a Hele-Shaw cell 

This flow is essentially two-dimensional, and one identifies the region occupied 
by fluid with a subset of the plane . In this paper we will consider the case of a 
viscous incompressible fluid that is surrounded by air while new fluid is injected 
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at a constant rate at the origin. If p is the pressure in the fluid then the velocity 
V of the fluid at that point (or really the mean velocity along the gap between the 
plates) is calculated as 

V = -Vp. (1) 

Here we have neglected some physical constants. The pressure in the fluid is har- 
monic except at the origin, where it will have a logarithmic singularity, and is zero 
on the boundary since we assume that both the air pressure and surface tension is 
zero. 

To model this, let {fit} be a family of domains in containing the origin, 
and let pt denote the function which is zero on dQt and which solves 

Apt = So 

on fit , where So is the Dirac measure at the origin. Then Q,t is called a solution to 
the (classical) Hele-Shaw flow if for all t 

Vt = -Vpt ondQt, (2) 

where Vt is the normal outward velocity of dilf Clearly we need some regularity 
of the family {f2t} in order for this to make sense, thus the question of local exis- 
tence of a classical solution to the Hele-Shaw flow with a given initial domain fio 
is non-trivial. 




Figure 2: Dynamics of the Hele-Shaw flow 

Similarly one can model the case of a Hele-Shaw cell where the fluid moves 
through a porous material. If the permeability of this material is given by a posi- 
tive function k then, by Darcy's law, the equation of motion of a fluid in the cell 
becomes 

Vt = — AcVpt outfit, (3) 

with Pt defined as before. In 1251 this situation is considered, where it goes under 
the name of elliptic growth of Beltrami type. It is also equivalent to studying 
the classical Hele-Shaw flow on a Riemann surface, with the Riemannian metric 
encoding the permeability (see f22|). 
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There is a vast literature on the Hele Shaw flow (see f20l and the references 
therein). The main short-time existence and uniqueness result for classical solu- 
tions to the Hele-Shaw flow in the case k = 1 is due to Kufarev and Vinogradov 
f401 and dates from 1948. It states that if flo is a simply connected domain with 
real analytic boundary, then there exists a unique (indeed real analytic) solution 
{Qt}t£{-6,e) to the Hele-Shaw flow for some e > 0. Observe here that the solu- 
tion extends both forward and backward in time. In |30| Reissig and von Wolfers- 
dorf gave a new proof of this result using a non-linear version of the Cauchy- 
Kovalevskaya theorem due to Nishida. Tian |38| provided yet another proof re- 
lying on properties of the Cauchy integral of the free boundary, and recently Lin 
1281 proved the same result using a result of Gustafsson on rational solutions 1181 
combined with a new perturbation theorem to get to the the general case. 

There is also a slightly different setting of the Hele-Shaw flow, where the 
source of fluid is not a point but a curve inside the starting domain (see e.g. 
l2llI0|[nj ). In LL3I Escher and Simonett proved short-time existence of classi- 
cal solutions to the Hele-Shaw flow in this setting, and in II II II2I they proved 
short-time existence of classical solutions to the Hele-Shaw flow with surface ten- 
sion. There is also related work of Hanzawa |2I| on classical solutions to the 
Stefan problem. 

Less has been written about the case of varying permeability. In |22| Heden- 
malm-Shimorin proved short-time existence under the assumption that n is real 
analytic, and the starting domain has real analytic boundary. They also proved 
short-time existence when the starting domain is empty, under an additional as- 
sumption that A log K > 0, which translates to saying that the Riemann surface in 
question has negative curvature. See also [23 1 for related results. 

2 Summary of Results 

Our first result is the short term existence and uniqueness of the Hele-Shaw flow: 

Theorem 2.1. Let Qo be smooth Jordan domain containing the origin and let k be 
a smooth positive function defined in a neighbourhood of dO,o- Then there exists 
a smooth positive function k' defined in a neighbourhood of dQ,a with k = n' 
outside of Qo , and an e > such that there exists a smooth increasing family of 
domains {Q,t}ts{-e,e) which solves the Hele-Shaw fiow with varying permeability 
k' . The extension k' is not unique, but a given choice uniquely determines the fiow, 
in the sense that any two such flows agree on their common domain. In particular 
we have short-time existence and uniqueness of the Hele-Shaw flow with varying 
permeability forward in time. 

We remark that this goes beyond the results of Hedenmalm-Shimorin since 
we do not require k or the domain of the initial domain to be real analytic. In 
fact if K = 1 then it is known the boundary of the domain being real analytic is 
necessary for there to be a solution backwards in time (see |39 | or 1 19, Thm 10]). 
Thus part of the above theorem says that this can be circumvented if we modify 
the permeability appropriately. Moreover if we do make this assumption we can 
say more and in particular recover Hedenmalm-Shimorin's existence result: 

Theorem 2.2. Suppose in addition to the previous assumptions that k is real an- 
alytic and f2o has real analytic boundary. Then we can take k — k in a neigh- 
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bourhood of dQo: ond so we have short-time existence and uniqueness of the Hele- 
Shaw flow both forward and backward in time. 

In our opinion, the main interest of this work is not the particulars of the above 
theorems but the techniques introduced for their proof. We shall show that a solu- 
tion to the Hele-Shaw flow is equivalent to being able to lift fit C = C to a 
family of holomorphic discs St in C x such that 

1. f2t is the image of the projection of St to C and 

2. The discs St attach along their boundaries to a certain totally real submani- 
foldAofC^ C C X P\ and 

3. The closure of St intersects C x {oo} C C x only at the point (0, cxa). 




Figure 3: Lifting the Hele-Shaw flow 



This correspondence comes about through thinking of the lifts St as the graph 
of certain "Schwarz" functions on Sit that are holomorphic except for a simple 
pole at the origin. Starting with t = 0, we construct a smooth strictly subharmonic 
function (j) (that encodes the permeability by k = 1/ A</>) such that f2o admits such 
a Schwarz function. Then setting 

A = {(.,f^)}cc^ 

oz 

the graph of So will be a holomorphic disc So whose boundary lies in A. 

Once this is done we are in a position to apply the well-developed theory of 
embedded holomorphic discs to prove that the moduli space of embedded discs 
nearby So that are attached to A and that pass through (0, oo) is a smooth manifold 
of real dimension one which gives both the existence and uniqueness results we 
require. This dimension count in turn relies on a standard argument using the 
associated Schottky double and a simple intersection count. 

Along the way we use the same technique to analyse the moduli space of 
quadrature domains (which will be defined in Section[5|. 
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Theorem 2.3. The moduli space of smooth quadrature domains of order n and 
connectivity c is a smooth manifold of real dimension 4n + c — 2. 

It is also interesting to consider the Hele-Shaw flow with empty initial con- 
dition. Despite this having a rather natural physical interpretation, there appear 
to be rather few results in this direction. This may be because many of the tech- 
niques in the study of the Hele-Shaw flow rely on identifying the starting domain 
with the unit disc (say through the Riemann Mapping Theorem) and thus cannot 
be applied. We shall show how the lifting technique above can be used to prove 
short term existence in this setting as well, which will require more sophisticated 
machinery concerning holomorphic discs, specifically results of Donaldson con- 
cerning a suitable homogeneous Monge Ampere equation. 

Theorem 2.4. Let k he a positive smooth function. Then there is an e > and 
family fit of domains in C containing the origin that satisfies the Hele-Shaw flow 
with permeability k and such that 



(i.e. the limit ofO,t as t tends to zero is just the origin). 

We remark that this improves on the results of Hedenmalm-Shimorin in that 
we do not need to assume k is real analytic or that A (log k) is positive. 

This lifting interpretation of the Hele-Shaw flow through Schwarz functions 
comes about through looking at the complex moments 



where dA denotes the Lebesgue measure on C. It is known that the Hele-Shaw 
flow is characterised by the condition that Mk {t) is constant with t for fc > 1 (and 
so the area Afo(t) is linear with t). One can similarly study the problem of flows 
Q,t for which Mk (t) for fc > 1 are allowed to vary in a particular way, and this puts 
the Hele-Shaw flow into the more general picture of "inverse potential problems". 
We use the same approach as above to prove a statement about these more general 
flows: 

Theorem 2.5. Let Q be a smooth Jordan domain, and V be a smooth, nowhere 
vanishing outward pointing normal vector field on Then there exists a smooth 
variation Sit for t £ [0, e) for some e > such that 



and whose moments vary linearly in t. 

Note that this theorem is not immediately obvious since a-priori there could 
be relations among the higher moments; thus we think of it as showing a kind of 
independence among them. In fact it is known that if k is real analytic, and Q, 
has real analytic boundary, then the moments Mk for fc > provide parameters 
for the space of domains with analytic boundary near to Q.. We refer the reader to 
Theorem |7.2| for a more precise version of Theorem |2.5| in this setting. If we set 
V — — Vp in the above theorem, where p denotes the Green function of SI, we 
recover the Hele-Shaw flow, and thus this result generalises Theorem |2.1| 



Area(rit) =t forO <t < € 
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Terminology: We identify with the complex plane C in the standard way. A 
domain C is a set that is open and connected and its boundary is dfl — Q — Q, 
where the bar denotes the topological closure. The connectivity of fl is the number 
of connected components of dfl. If dO, can be written locally as the graph of a 
smooth function then we shall say that Q has smooth boundary or simply that SI 
is smooth, with analogous definitions if this graph can be taken to be real analytic, 
and in this case there is a well defined unit outward normal vector field n on d^l. 
We say S7 is a smooth Jordan domain if it is the interior region determined by a 
smooth Jordan curve (i.e. a smooth non-self intersecting loop in R^). 

If 7 C R is an interval then a family of domains fit for t £ I is increasing if 
fit C flti for t < t' and the family fit is said to be smooth if the boundary dQt can 
be written locally as the graph of a smooth function that depends smoothly on t. 
If flo is smooth with unit normal vector field n on OQq then an increasing smooth 
family Qt for t £ (— e, e) is determined by dQt — {x + fix, t)nx : x G dQ.o} for 
some positive smooth function ,ft{x) = f{x, t) on 9f2o, and the normal outward 
velocity of fit at t = is 

3 Lifting the Hele Shaw Flow 
3.1 Complex moments 

Let fit be a family of domains in R^ = C and k be a positive smooth function. 
For fc > the fcth complex moment of fit is 

Mfc(t) ;= / z''-dA, 
Jilt 1^ 

where dA denotes Lebesgue measure. Clearly then Mo{t) is just the area of fit, 
and we refer to the Mk for fc > 1 as the higher moments. 

An important discovery by Richardson f31| is that these moments are con- 
served by the Hele-Shaw flow. Since we are assuming the fluid is injected at 
a constant rate, we may as well reparametrise t so that Mo{t) — Area(fif) — 
Mo(0) + t. The upshot then is that this system is integrable, by which we mean 
that these moments characterise the Hele-Shaw flow: 

Theorem 3.1. 

1. Let fit be given by the Hele-Shaw flow. Then the higher moments Mk (t) 
remain constant with t. 

2. Conversely, any smooth increasing family of simply connected domains { fit } 
for t £ (— e, e) that contain the origin and whose higher moments Mk are 
constant with respect to t is a solution to the Hele-Shaw flow ( after a possible 
reparameterising oft). 
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The first statement is Richardson's calculation, namely that if the Hele-Shaw 
equation (|3} is satisfied then 

J sit J dUt 

which vanishes if fc > and equals 1 when A: = as pt = on d^lt and A{pt) — 
—Sq. The second statement says in effect that this calculation can be reversed (see 
l20l p22] or L22J ). 

3.2 Schwarz functions 

We next discuss how the moment characterisation of the Hele-Shaw flow can be 
recast in terms of the existence of Schwarz functions. To describe this, fix a pair 
(51, (f)) where 57 is a domain in — C, and </) is a smooth function defined in some 
neighbourhood of d^l that is strictly subharmonic (i.e. such that A<j) is positive). 

Definition 3.2. We say that a function 5 is a Schwarz function of (57, (j>) if S is a 
meromorphic function on 57 which extends continuously to 957 and such that 

S{z)^^^ on 957. 

Remark 3.3. The model case usually considered is where (f> — \z\^, so S{z) — z 
on 957. We remark that what we define above is sometimes referred to as an interior 
Schwarz function, since we are requiring that it be defined on all of 57. Much of 
the work of Schwarz functions asks instead that S be defined, and holomorphic, 
on some neighbourhood of 957. It turns out, in the model case 4>{z) = \z\'^, that 
such an 5* exists in this restricted sense if and only if 57 has real analytic boundary 
(see e.g. the book 1 37J by Shapiro), and extends to a meromorphic function on 57 
if and only if 57 is a quadrature domain jT). 

The connection between Schwarz functions in the model case (f>{z) = l^l^ and 
the Hele-Shaw flow is well known (e.g. L20. 3.7.1]). The following generalises this 
to suit our needs: 

Proposition 3.4. Let (j) be a smooth and strictly subharmonic function defined in 
a domain D. Suppose 57f, t G (fl,&), is a smooth increasing family of Jordan 
domains containing such that for all t £ {a, b), 957t C D, and for each t there 
exists a Schwarz function St of the pair (57t, that is holomorphic on 57t except 
for having a simple pole at 0. Then after some reparametrization s{t), the 57s are 
a solution to the Hele-Shaw flow with varying permeability k = 1/ Atfi. 
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Proof. Pick to G (a, b) and let t £ {to,b). Then for fc > we have, 



/ 



if kd(p if kd(j) 

z -—dz / z -^dz 



2 J dz 2 ign^^ dz 

if k r< 1 if k n 1 

= - / z btdz— - / z btgdz 
^ J ant ^ Jaritg 

= ^ z^BStdz-^ I z'^dStgdz^ 

= z''(0)nReso{St~ Sta). 

Thus the complex moments are preserved, and so fit is a classical solution to the 
Hele-Shaw flow ( |3.1[ l (perhaps after reparametrization). □ 

Thus we wish to prove the existence of Schwarz functions St, which is made 
easier through the freedom of choice in the subharmonic function (fi. Our next 
aim is to construct a suitable that solves this problem for t = 0. To guide 
the subsequent discussion we first recall a classical argument in the case where 
4> ~ \z\^ . Consider the Cauchy integral 



Define /+ as the restriction of / to Q,o and /_ as the restriction of / to the comple- 
ment of fio- Clearly /+ and /_ are both holomorphic in their respective domain of 
definition. A careful analysis shows that both /+ and /_ extend smoothly to d^lo 
and the Sokhotski-Plemelj formula states that on the boundary 

f+{z)^z + f-(z) (4) 

(see, for example, the book (151 by Gakhov). Moreover it can be shown that the 
boundary of Qq is real analytic if and only if both /+ and /_ extend holomorphi- 
cally to a neighbourhood of d^lo, and in this case S :— f+ — f- is holomorphic 
on this neighbourhood and S(z) — z on dO,. 

We state first a simple version of what we shall prove: 

Proposition 3.5. Let Qo be a smooth Jordan domain containing 0. Then there 
exists a smooth strictly subharmonic function (f) defined on a neighbourhood U of 
the complement ofQo such that A(j> = 1 on (7 — f2o cind so that (ilo, <t>) admits a 
Schwarz function So that is holomorphic except for having a simple pole at 0. 

For the problem of varying permeability we need something a little stronger: 

Proposition 3.6. Let Q,o be a smooth Jordan domain containing and n be a 
smooth positive function defined on a neighbourhood of d0.o- Then there exists a 
smooth function k defined in on a neighbourhood U ofdilo such that k' = k on 
U — ilo dnd a smooth strictly subharmonic function (j) defined on U with Acf> — 
1/k' such that (D,o,cf>) admits a Schwarz function So that is holomorphic except 
for having a simple pole at 0. 

Moreover if k is real analytic, and Q,o has real analytic boundary then we can 
take k! = K. 
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Proof. We need only prove Proposition ! 3. 6 1 assume k is defined 

on all of C and is smooth and positive there. Pick some smooth function (I>q on C 
such that A0O = 

Consider the Cauchy integral 



J ana 



900 



c- 



As above we write /+ and /_ for the restriction of / to f^o and to the complement 
of Slo respectively. Then the Sokhotski-Plemelj formula says that on 9Slo, 

/4.) = ^ +/_(.). (5) 

Observe next that /_ tends to zero as z tends to infinity. Thus for some constant 
a we have that 

= f +/», 

where / vanishes to second order at infinity. Such a holomorphic function pos- 
sesses a i9/9z-primitive F on C — S7o. Letting h = 2Re(F) we can rewrite 
equation 151 as 

f+{z) - z ^ Q^{4>o + h) on dflo. (6) 



z oz 



Now set 



4> ~ 'f'o + h on C — fio , 

and observe that since h is the real part of a holomorphic function, A(j> — A(j)o — 
1/ hi. Moreover since /_ had a smooth extension to dQ.o it follows that cf) extends 
smoothly to a neighbourhood of the complement of fio and we can define k' = 
l/{A(f)) which will be positive near dQ.o. Finally 5*0 := f+{z) — j is a Schwarz 
function of (fio, <i>), which is holomorphic except for a simple pole at 0. 

If K, is real analytic and dQ.o has real analytic boundary, then /_ extends to 
a holomorphic function over a neighbourhood U of dQ,o 1241 1.6.4,2.5.1], and h 
similarly extends to a harmonic function over the boundary. Thus defining <j) over 
this neighbourhood as above we have A0 = 1/ n over this neighbourhood as well, 
and so n' — k. □ 

Remark 3.7. From the above proof one sees that in fact the second statement 
holds under the condition that /_ extends analytically over dQ, the arguments 
below show that in this case there is a backwards solution to the Hele-Shaw flow. 

Remark 3.8. In what follows we will use in an essential way that (j) is smooth on 
a neighbourhood of d^lo- But one should note that the Sokhotski-Plemelj formula 
holds for much more general domains (e.g. Lipschitz domains) and one still gets a 
Schwarz function of (JIq , , in this case and i 



3.3 Holomorphic Discs 

We are now ready to phrase the Hele-Shaw problem in terms of holomorphic 
discs. Let <j) be the subharmonic function defined in Proposition |3.6[ and So be 
the Schwarz function for (f2o, 4>). Define 

A~{{z,^)}cCxC. 
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Since this is the graph of a real strictly subharmonic function, it enjoys various 
properties that will be discussed in the next section (in fact it is a totally real LS- 
manifold). We shall consider a meromorphic function on an open set in C as a 
function to the Riemann sphere = CU{oo}in the standard way, giving it the 
value oo at its poles. 

Deflnition 3.9. Suppose (SI, 0) admits a Schwarz function S holomorphic except 
for having a simple pole at 0. Define 

E = En~{{z,s{z))eCxr'- -.zen}. (7) 

So En is a holomorphic disc in C x . Moreover, 

1. If TTi denotes the projection to the first factor then tti : E — )■ vri(E) is an 
isomorphism, 

2. The boundary of E lies in A, 

3. The closure of E intersects {0} x P^ precisely in the point (0, oo). 

Theorem 3.10. Suppose that Et is a smooth family of holomorphic discs in C X P^ 
that satisfy the three conditions above and such that Qt ■= 7i"i(Ei) is an increasing 
family of smooth Jordan domains. Then Qt is a solution to the Hele-Shaw flow with 
permeability K := l/(A</i). Moreover any smooth increasing family Q,t of smooth 
Jordan domains that solves the Hele-Shaw flow with permeability k arises in this 
way. 

Proof. The first statement just summarises what we have said thus far. Given 
Et, it can be realised as the graph of a Schwartz function 5*4 for Q.t ~ 7ri(Et). 
Condition (3) says that St is holomorphic except for having a simple pole at 0. 
Thus the result we want follows from Proposition |3.4| 

For the second statement, suppose that Q.t is a smooth increasing family of 
smooth Jordan domains that satisfy the Hele-Shaw flow. Let pt{z) = p{z, t) be 
the pressure which is required to satisfies Apt = —So on Q,t and pt = on dQ.t. 
Now define 

ut{z) = u{z, t) 

and set 

^ _ d(f> dut 
oz oz 

where denotes the characteristic function of fio- Then = 1 — x^o ^ ^^o 
in the sense of distributions, (see [20i Sec 3.1] or Remark [3 . 1 2| below) and one 
easily checks that St is a Schwarz function for fit which is holomorphic except for 
a simple pole at . Thus the Q.t are induced from the corresponding family En^ 
as required. □ 

In the following section we show that the moduli space of such holomorphic 
discs has real dimension 1, and thus our original disc Eo deforms in a smooth one 
dimensional family Et. Observe that condition (1) and the property that tti (Et) is 
a smooth Jordan domain are open conditions, and thus since they hold for Eo they 
will hold for Et for small t as well. 



= / p{z,T)dT 



XfJo ( ^ - So 
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Remark 3.11. In much the same way one can interpret the Hele-Shaw tlow with 
multiple inputs, in which case condition (3) above must be suitable modified. We 
leave the details to the reader. 

Remark 3.12. To understand the long-time behaviour of the Hele-Shaw flow, and 
also for more general starting domains, one has a weak formulation of the problem. 
Indeed, following Elliott- Janov sky 1 10|, Gustafsson 1 19| and Sakai 1321 . 1351 it is 
possible to cast a weak formulation of the Hele-Shaw flow as an obstacle problem. 

Let / be a distribution on some domain D. A function u on D is said to solve 
the corresponding obstacle problem if 

- Au - / > 0, u>0, (-Alt - f)u = 0. (8) 

The free boundary T{u) of the obstacle problem is the boundary of the coincidence 
A(u) :— {x : u{x) = 0}. If the domain D has a boundary one also specifies a 
non-negative function g on dD and demand that the solution u should, in addition 
to satisfying ([8|, be equal to g on the boundary. 
Now let D = and for any t > set 

ft ■■= Xno -l + tSo. 

By general potential theory the corresponding obstacle problem has an unique 
lower semicontinuous solution ut and one calls the family 

fit — {z : ut{z) > 0} 

the weak solution to the Hele-Shaw flow with initial domain Slo • It is know that a 
smooth solution, if it exists, will coincide with the weak solution, and if the weak 
solution is regular enough, it is the smooth solution (see e.g. f20|). Moreover there 
does exists a weak solution that is defined for all time, but in general the f2t will 
eventually be non-smooth or cease to be simply connected. 

Now suppose 5*0 is a Schwarz function of {yio, 4>) and suppose that A(j!) — 1 
on the complement of Qq . We have seen that such a function exists at least when 
r^o is Lipschitz. Let Qt be the weak solution to the Hele-Shaw flow. We observe 
that for any t £ [0, oo), there exists a Schwarz function for (Qt, (j)), holomorphic 
except for a simple pole at 0. To see this, let Ut be solution to the above obstacle 
problem described. Then one simply notes that the function 

.dcf) dut. ,d(f> , 

St ■■= {-r. ~ XsioK-r. Oo) 

oz oz oz 

is the claimed Schwarz function. 

Hence using the same lifting we thus see that there is a natural deformation of 
our initial holomorphic disc Eo that exists for all positive time, but whose members 
need not be smooth and may have more exotic topology. However the possibilities 
are constrained, as there are strong regularity results due to Sakai I33i i34i for the 
boundaries dQ.t of the domains Sit belonging to a weak Hele-Shaw flow which 
says that if Q.(, C f2j (which, for example is the case if dQ,o is C^) then dQ.t is 
real analytic except for having at most a finite number of cusps or double points. 
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4 Holomorphic curves with boundaries in a to- 
tally real submanifold 

We have thus seen that the Hele-Shaw flow can be understood from the point of 
view of Schwarz functions, and in turn in terms of holomorphic discs with bound- 
aries in a given totally real submanifold. There is a rich literature on this topic, we 
now summarise parts that are relevant to our situation. 

4.1 The moduli space 

Let Af be a complex manifold of complex dimension n. A submanifold A is said 
to be totally real if TA n J{TA) = {0}, where J denotes the complex structure 
on M. A totally real submanifold is said to be maximal if it has real dimension n. 

Let (E, 3E) be a Riemann surface with boundary. A holomorphic non-const- 
ant map / : E — >■ M is called a holomorphic curve in M if / extends to a 
function on 9E. If / is an embedding then we call the image of / an embedded 
holomorphic curve. By abuse of terminology we shall refer to this image as E, and 
say that it has boundary in A if /(9E) C A. In the special case that E is the unit 
disc D C C and f{dD) C A we call the image an embedded holomorphic disc 
with boundary in A. By the smooth Riemann mapping theorem |7 Thm 3.4], one 
can equivalently replace the disc D with any smooth Jordan domain. 

The moduli space of holomorphic curves with boundary in a totally real max- 
imal submanifold has been studied by numerous authors in different contexts, e.g. 
Bishop 1 5 1, Donaldson |9|, Forstneric 1 14] and Gromov 1 16]. In complex analysis 
this is related to polynomial hulls and the homogeneous Monge-Ampere equation 
L9J , while in symplectic topology it is used in the study of Lagrangian submani- 
folds II6I . More recently LeBrun and Mason 1261 1271 have shown a connection to 
twistor theory. 

The first step in this theory is to give the space of all (not necessarily holomor- 
phic) embedded curves in some regularity class with boundary in A the structure 
of a Banach manifold B. Let E be such an embedded curve, let TV denote the nor- 
mal bundle and let E be the normal bundle of 9E in A. We can identify E with a 
totally real subbundle of A'^ restricted to 9E. Locally around E the Banach mani- 
fold B is modelled on the Banach space C'°+^'"(E, N, E) of sections of N with 
boundary values in E (this construction uses the exponential map with respect to 
some appropriately chosen metric, see |27|). 

The next step is to consider a Banach vector bundle n : £ ^ B whose fibre 
over E is 

C'='"(E,A'''' ®Ar) 

and a canonical section S : B ^ £ whose zero locus M := 5"^(0) is the 
space of holomorphic curves. The implicit function theorem for Banach spaces 
says that is a smooth manifold provided that the graph of S is transverse to 
the zero section. This is so if and only if at any point E in A^, the linearisation 
-De, which by definition is DS composed with the projection to the fibre of £, 
is surjective and has a right inverse. Now De is nothing but the Cauchy-Riemann 
operator B that maps sections of A'^ with boundary values in _E to (0, 1) -forms with 
values in A'^. It is well-known that this operator is Fredholm, so what remains to 
establish transversality of S is surjectivity. Thus if the Cauchy-Riemann operator 
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is surjective, then Mlisa. smooth manifold and its dimension is given by the index 

of a. 

In general the problem of proving transversality involves picking a new generic 
almost complex structure on the manifold M. However for our purposes M has 
complex dimension 2 and then the problem becomes significantly easier and, as we 
shall see, reduces to the computation of a single topological invariant of a related 
holomorphic curve. 

4.2 Dimension count and the Schottky double 

Following LeBrun in 1271 . we now discuss how to use the Schottky double to 
compute the dimension of the moduli space A4 of embedded holomorphic curves 
in a manifold A/ of complex dimension 2. 

Let E be a compact Riemann surface with boundary SE. By gluing it together 
with its complex conjugate S along SE we get a closed Riemann surface X (as- 
suming some regularity of E, see 1271 '). In the theory of quadrature domains E is 
called the Schottky double of E fTTj- The complex conjugate iV of A'^ is a holo- 
morphic vector bundle on E, and by gluing and A'^ along dT, so that E = E we 
get a holomorphic vector bundle on X. 



There is an antiholomorphic involution pon X which lifts to A/" such that E is 
the fixed point set of p. As noted in li27J . this implies that E is real analytic as a 
subspace of M. This involution thus acts on each of the vector spaces 



C'=+'^°(A,0(AA)) and H''{X,0(U)) forp = 0, 1. 
Let Cp+^-"{X, 0(J\f)) and H^iX, 0{J\f)) denote the invariant subspaces. Thus 
H'\X, 0{N)) = Hl{X, 0{N)) e iHliX, 0{N)). 
LeBrun shows in 1271 that 



N 




X N 



Figure 4: The Schottky double 



(E,iV, E)^c; 



,k + \,a 



{X,0{N)) 



with a canonical isomorphism. Also f27' Lemma 1] the kernel of the Cauchy- 
Riemann operator is canonically isomorphic to Hp(X, 0{Af)) while the cokernel 
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is canonically isomorphic to Hp{X ,0{Af)). Thus d is surjective if and only if 
h^{X, 0{Af)) = 0, and then the index is equal to h°{X, 0{J^)). 

Now suppose that when the complex dimension of the ambient manifold M 
is two. Then A/" is a complex line bundle, and so as long as degA/" > 2g — 2, 
where g is the genus of X, then the numbers h° {X, 0{N)) and {X, 0{Af)) are 
topological invariants that only depend on g and the degree of Af. 

If A happens to be the fixed point set of an antiholomorphic involution p of 
M, then the Schottky double A' of E can be identified with the Riemann surface 
S U and A/" is then identified with the normal bundle of S U in AI. 

Hence a strategy to compute the numbers h^\X, 0{Af)) and {X, 0{N)) for a 
general pair (S, A) in the four dimensional case is to continuously perturb (E, A) 
to some (E', A') where A' is the fixed point set of an involution p, and then use the 
adjunction formulae to compute the degree of the normal bundle of E' U p(E'). 

5 Moduli spaces of quadrature domains 

We digress in this section from the Hele-Shaw flow, and consider the above setup 
in the case A = {[z,z)} which is related the concept of quadrature domains. 

Deflnition 5.1. A bounded domain f2 in C is called a quadrature domain if there 
exist finitely many points ai, am £ Q. and coefficients Ckj such that for any 
integrable holomorphic function f onQ, the quadrature identity holds: 

fc = l 3=0 

The integer n = ^JIL-^ is called the order of the quadrature domain (assuming 
Cfc,nfc-i 7^ 0). We say that a quadrature domain is smooth if its boundary is 
smooth. We also recall that the connectivity c of a domain is the number of 
connected components of dO,. 

Assume that [Q,, 4'{z) — \z\'^) has a Schwarz function 5*. Then 

/ fdA^'- [ fd{-zdz) = I f fzdz 
Jn ^ Jn ^ J an 

= ^ [ fSdz [ fdSdz 
2 Jan 2 Jq 

m — 1 
fc = l j=0 

where at are the poles of S and nt is the order of the pole at at, and thus is a 
quadrature domain. It is proved in (Tl Lemma 2.3] that the converse holds as well. 
Clearly then the order of the domain n equals the number of poles of 5* counted 
with multiplicity. 

It is also clear that Schwarz function S is uniquely determined by the domain, 
since if we had two Schwarz functions, their difference would be a meromorphic 
function that was zero on the boundary, and therefore this difference must be zero. 

From now on assume that the quadrature domain Q, has smooth boundary. Con- 
sider the graph 

E — {iz,s{z)) -.zen} 



SHORT-TIME EXISTENCE OF THE HELE-SHAW FLOW 



as an embedded holomorphic curve in C x . It is attached along its boundary to 
the totally real submanifold 

A = {{z,z) -.zeC}. 

Since A is the fixed point set of the antiholomoiphic involution p : (z, w) i— >■ 
(«), z) , the Schottky double X is naturally identified with S U p(E) and the normal 
bundle A/" is identified with the normal bundle of E U p(E). 

We wish to calculate the degree of Af. By the adjunction formula, 

where 0{X) denotes the line bundle associated to the divisor X. The line bundle 
Kfi is represented by the divisor —2{H\ + H2) where H\ and H2 denote the 
hyperplanes x {00} and {00} x P^ respectively. Thus 

degAA = 2{Hi + H2)-X + degKx. 

where the dot represents the intersection number of the relevant divisors. Clearly 

Hi ■ X ^ H2 ■ X ^ n, 

and 

degKx ^2g- 2, 
where g denotes the genus of X. Combined this yields that 

degA/" = 4n + 2g - 2. 

Now n > 0, so degA/" > 2g + 2 which implies h^{X,J\f) = 0, and by Riemann- 
Roch that 

h^iX,^/) = degAA-g + 1 = 4n + .g- 1. 
Finally note that the genus g is precisely c — 1, where c is the connectivity number 

of n. 

Hence the moduli space of embedded holomorphic curves with boundaries in A 
is a smooth manifold of real dimension 4n + c — 2 near E. Any such holomorphic 
curve lying close enough to E gives a Schwarz function of a smoothly bounded 
quadrature domain. By the fact that all quadrature domains have unique Schwarz 
functions (again with respect to (j>{z) — \z\'^), it follows that the moduli space 
of smoothly bounded quadrature domains of order n > and connectivity <? is a 
manifold of real dimension 4n + c — 2 . We have thus proved Theorem |2.3| as stated 
in the introduction. 



6 Short-time existence of the Hele-Shaw flow 

This section completes the proof of Theorem |2.1| To ease notation let = ilo 
be a smooth Jordan domain. We saw above that we could construct a generalised 
Schwarz function S of the pair {Q, (j)) where A<j) = 1/k on the complement of Q 
and (j) extends smoothly to a neighbourhood of the complement. If E denotes the 
graph of S and A the graph of ^ we have that E is a holomorphic disc whose 
boundary lies in the totally real submanifold A. 
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We want to analyse the moduli space of holomorphic discs near E with bound- 
ary in A that intersects the hyperplane Hi = x {00} at the point (0, 00). 

First we consider the moduli space of all such discs, irrespective of their in- 
tersection with Hi. We need to calculate the numbers h°{X,jV) and h^{X,jV). 
As we continuously deform the domain fl to the unit disc Q', the pair (A", A/") and 
submanifold A continuously deform to the pair {X' ,Af') and to A' — {(2, z) : 
z £ C} respectively. Thus the Schwarz function for fi' is simply -^j i.e. the usual 
(interior) Schwarz function of the unit disc thought of as a quadrature domain. 
Since the numbers h^\X,Af) and h^{X,Af) are invariant under continuous de- 
formation, the calculation in the previous section yields that h^{X,J^) = and 
h^\X,jV) = 3, since n — 1 and g = 0. Thus the required moduli space is a 
smooth manifold of real dimension 3 near E. 

Considering moduli spaces of holomorphic curves passing through a given 
point is also standard in the theory (see, for example, 1 29 1). We simply replace our 
Banach manifold B with the subspace of curves going through the point, which 
is locally modelled on the Banach subspace of sections to the normal bundle that 
vanish at that point. We will denote this space by C'^"*"^'"(E, A*', E)o. Following 
1271 it follows that this space is naturally isomorphic to the space of p invariant 
sections of A/" that vanish at both (0, cxd) and p((0, 00)). We denote this space by 
Cp^^'" {X ,jV)o- The section 5 and its linearisation 8 are simply the restrictions 
of the original ones. 

We check first that d is still surjective. Pick an arbitrary element 

g G C'='"(E,A"'i ® Af). 

By the above there exists a section / £ C'°+^'°'(E, TV, _B) such that Bf = g. 
Since the degree of Af is positive and X is isomorphic to we know that there 
is a holomorphic section, s say, of A/" that is equal to / at (0, 00) but zero at 
p((0, 00)). Then the section s + p(s) lies in Hp{X, N), and thus corresponds to 
an element s in the kernel of d. We have that s = / at the point (0, 00) . Thus / — s 
lies in C'°+^'"(E, TV, E)o and d{f - s) ^ g. This shows that the restriction of 8 
is still surjective. 

The kernel of the restricted Cauchy-Riemann operator is of course the intersec- 
tion of the original kernel with the subspace C*^^'"(E, TV, E)o. Under the natu- 
ral isomorphism this equals the intersection of Hp(X, TV) and Cp^^'" (A", A/')o. 
Let H^{X, TV)o denote the subspace of holomorphic sections that vanish at both 
(0, 00) and p((0, 00)). The involution p acts on this space by 

p{s){x) ■- p{s{p{x))), 

and if Hp{X, TV)o denotes the invariant real subspace we have that 

H'\X, N)o = H^piX, N)o e iH°p{X, N)o. (9) 

Clearly_/f°(A',A')o is the intersection of H°{X,N) and Cp+^'^iX ,jV)o, and 
from |9|l we get that the real dimension of Hp{X, TV)o is equal to the complex 
dimension of H°{X,N)o. Since the degree of A/" is two and X the complex 
dimension ofH°{X,N)o is one. Thus the index of the restricted Cauchy-Riemann 
operator is one. 

Together with the surjectivity this implies that the moduli space of nearby holo- 
morphic discs with boundary in A and going through (0, 00) is a smooth manifold 
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of real dimension one. Thus Eo = S varies in a smooth family St for t suffi- 
ciently small, and by Theorem |3.10| it remains to prove that the family of domains 
r^t = 7r(Et) is increasing. 

The tangent space of the family of discs at E is given by Hp{X, N)o restricted 
to E. Since a non-zero element in Hp{X, N)o vanishes only at the points (0, oo) 
and p((0, oo)) it is nonvanishing along the boundary of E. The real subspace E 
is one-dimensional and is has one direction going out of Q,o and one direction 
going in. The nonvanishing implies that a nonzero tangent points either outwards 
everywhere or inwards everywhere. Thus (replacing t with —t if necessary) the 
family fit is increasing. 



7 The inverse potential problem and moment 
flows 

Let p be a density function on C (so we have changed notation from previous 
sections in which this density was played by 1/ k). The associated potential of a 
domain is the subharmonic function 



U^{z) := / log \z — w\'^ p{'w)dA{'w) 
Jn 



where dA is the Lebesgue measure. Note that is harmonic on the complement 
of Q,. Let D be some big disc containing Q,. The inverse potential problem asks 
whether a given harmonic function U on the complement of D (with proper growth 
at infinity) can be realised (uniquely) as the potential of a domain flin D. 

To see the relationship with the complex moments consider 

G:= — , 

which is holomorphic if U is harmonic. If f/ = for some domain O, then 

G{z) = -^ [ log\z-wfp{w)dA{w)= f ^^dA{w). (10) 
oz Jq Jq z — w 

Now for large |z| we can write this integral as a convergent power series m z~^ , 



where 



/ P^dA{v:)^y^a,z-\ 



flfc := I w'' ^ p{'w)dA{w) 
n 



Thus the coefficients at of the power series of G at infinity equals the complex 
moments of Q with respect to the measure pdA. Thus the inverse potential problem 
be restated as asking for a (unique) domain with prescribed complex moments. 

This potential problem has a rich history, both for the plane and for higher di- 
mensions, and also for more complicated kernels replacing the expression In | z — 
above (see, for example, |24| and the references therein). For this particular 
problem in the plane it is known that prescribing the complex moments does not 
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uniquely determine the domain (although it appears to be unknown if this unique- 
ness can fail for domains with smooth boundary). However starting with a suitable 
smooth domain Qo, for small variations Ut of the potential U^" there is a unique 
Qt close to Q such that Ut = U^* . Thus one can think of the set of moments as 
giving local parameters for the space of domains. We refer the reader to |24 Sec 
2.4] and |6] Sec. 2] for precise statements. We also point out the interesting work 
of Wiegmann and Zabrodin in this context, which connects these parameters to an 
integrable system coming from the dispersionless Toda lattice |41 1. 

We have seen that the Hele-Shaw flow can be phrased as a local existence 
question related to the inverse potential problem. Namely, if is a smooth Jordan 
domain containing the origin, then there is a non-trivial one parameter family of 
smooth domains Qt, all having the same higher moments as and whose area 
increases linearly. In this section we will prove short-time existence of a similar 
flow in which all the higher moments are allowed to vary linearly. Variants of the 
Hele-Shaw flow under some other force (e.g. gravity) are considered by Escher- 
Simonett 1111 1121 [T3J and this gives a physical interpretation to some of these 
variations. 

Theorem 7.1. Let Q be a smooth Jordan domain, and V be a smooth, nowhere 
vanishing outward pointing normal vector field on dil. Then there exists a smooth 
variation Q,tfort G [0, e) for some e > such that 



whose moments vary linearly in t. 

Note that without the assumption on the vector field V pointing outwards this 
theorem is not true, for example in the Hele-Shaw flow when p = 1 and V = Vp, 
then this would mean solving the Hele-Shaw flow backwards in time, which is 
impossible unless the boundary of Q. is real analytic. 

Of course the way in which the moments vary will depend on V . To discuss 
this suppose V is any normal vector field defined on the boundary of a smooth 
Jordan domain Q.. Write the measure V pda as the restriction of a complex one 
form g{z)dz to the curve dil. The complex-valued function g is then uniquely 
determined on d^l. Consider the Cauchy integral 



and as above let /+ = ' denote the restriction of / to and /_ = /_ ' 
the restriction to the complement of Q. Recall that both /+ and /_ has a smooth 
extension to the boundary and 



Now if Qt is a smooth variation of Q with initial direction V, the variation of 
the complex moments Mkit) of fit satisfy 




U=9 + f- 



on do.. 
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Since /+ is holomorphic in D, it then follows that 

dMk, 



dt 



|t=o = - z''f^{z)dz. 
J an 



Thus writing /_ as a power series in a neighbourhood of infinity, 

oo 
fc = l 

we conclude 

dMk. „ ., 

|f=o ~ —2-Kibk+i. 

We think of the condition that extend to an analytic function across dfl as a 
kind of "growth condition" on the coefficients bk that determine the variation of 
the moments. With this said we can prove something more precise: 

Theorem 7.2. Let Q he a smooth Jordan domain. Suppose complex coefficients 
bjc for k > are given with ibo G M, so that the function /_ = J]]^. bkZ~ 
defined in a neighbourhood of infinity extends to a holomorphic function across 
do,. Then there exists a smooth density function p' with p' = p outside of Q, and 
a smooth family {Q,t}tii(-t,i)< whose moments Mk{t) taken with respect to the 
measure p'dA vary linearly as 

Mk{t) = Mfe(O) - 2-Kihk+it. 

If moreover f2 has real analytic boundary and p is real analytic then we can take 
p'^p 

Remark 7.3. This theorem gives one precise way to state that these moments form 
parameters for the space of domains near SI. The final statement, under the extra 
assumption that the boundary is real analytic, appears to be well known and to have 
been proved in a number of contexts; it follows, for instance, from fS' Sec. 2]. We 
are not aware of a similar statement in the smooth case, or one such as Theorem 
|7.1| that uses the extra data that the initial normal vector field is outward pointing. 

The proof starts with the following uniqueness result: 

Proposition 7.4. Let Vi and V2 be two smooth normal vector fields defined on the 
boundary ofQ,. If for all k > we have that 



Ian 

then it follows that Vi = V2 



/ z^V\pds = / z'^V2pds, 
Jan Jan 



(11) 



Proof. Let gi and g2 be the complex valued functions so that Vipds = Qidz on 
9n for i = 1,2. By our calculations above it follows from \\\\ that gi — 32 
extends holomorphically to Q.. Let it be a conformal map from the unit disc to Q, 
and consider the pullback of (gi ~ g2)dz to the unit disc. It will be a holomorphic 
one form which can be written as h{w)dw, where w denotes the holomorphic 
coordinate on the unit disc. The fact that (31 — g2)dz restricts to a real measure 
on 9n implies that the restriction of h{w)dw to the unit circle is real. Thus in turn 
implies that wh{w) is real valued on the unit circle. It follows that 'wh{w) = 
and thus gi = 32 , and consequently Vi — V2. □ 
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Recall that we have a totally real submanifold A which is the graph of 
where — p outside of il. We also have a holomorphic disc S attached to A 
which projects down to As we saw in Section [6] E belongs to a one parameter 
family of holomorphic discs attached to A which goes through the point (0, oo), 
and the projection of these yields the Hele-Shaw flow (the permeability being 1/p). 

We record the following theorem. 

Theorem 7.5. Let E and A be as above. Given any small smooth perturbation As 
of K, the family Et deforms to a smooth family Tia,t of holomorphic discs going 
through (0, oo) attached to Aj, along their boimdaries. 

Proof. This is standard in the theory of moduli spaces of embedded holomorphic 
curves (see e.g. 1271 '). Let A^ be a smooth parametrised perturbation of A, s G 
(— £, e). The Banach space B one considers is the space of embedded discs in 
some regularity class together with the parameter s, where the discs passes through 
(0, oo) and is attached to A^. The Banach bundle £ is defined just as before, and 
so is the section S. The linearisation of S is the B operator which maps the tangent 
space of B to the fibre of £. Let Es be a smooth family of smooth discs lying in B 
such that Eo = E and Es attaches to A^. Then the tangent space of B is spanned 
by C'°+^'"(E, N, E)o together with the derivative of E^ at s = 0, which we will 
denote by /. We have already established that the linearisation d is surjective, even 
when restricted to the subspace C*"''^'"(E, N, E)q. Therefore there is a section g 
in C*+^'"(E, iV, E)o such that df = Bg. Thus the kernel of B is spanned by the 
kernel of the restricted B operator together with / — g, so it has dimension 2. By 
the implicit function theorem for Banach manifolds the space of holomorphic discs 
nearby E lying in B is a manifold of real dimension 2. Also the subsets which are 
attached to Aa for a given s will be one dimensional submanifolds. □ 

We now prove Theorem |7.l| Let /_ be the restriction to the complement of Cl 

f=f ^dy., 
Jan z-w 

where Vpds = g{z)dz, and choose some smooth extension to a neighbourhood U 
of the boundary of SI. We will consider the smooth perturbation of A given by 

A, ~{{z,'^ + sf.):z&UZ)Bn}. 
oz 

By Theorem |7.5| the family of holomorphic discs corresponding to the Hele- 
Shaw flow Ef perturbs to a smooth family E^.t with boundary in As. For each 
s select the holomorphic disc Eg := Es,t(s) so that the corresponding Schwartz 
function Sa has the same residue as So. That Es is a smooth family follows from 
the ordinary implicit function theorem in finite dimensions. Let ils be the projec- 
tion of Es, we have for small s, 

f z'^pdA— I z'' pdA = f z'^^^dz— f z^^^dz = 

z^{Sa- sf-)dz- [ z^Sodz = -s[ z'^f-dz. (12) 
Jan^ J asiQ Jans 

Here we used the fact that •S's = §| + sf- on dQ,s and that Ss and 5*0 have the 
same residue at 0. 
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If Vb is the initial normal vector field on 957 corresponding to the perturbation 
fis, then \12\ implies that Vo and V induce the same initial variation of moments. 
From Proposition |T4] it thus follows that the initial direction of Q.s is given by V. 
In particular this initial direction is, by assumption, non-vanishing and outward 
pointing, so by smoothness dQ.s will lie in the complement of Cl for small s. But 
on this complement /_ is holomorphic, and 

/ z'' f-dz = 2TTibk+i 
J aOs 

(which in particular is independent of s). Thus by ^12\ again, all the moments of 
fla vary linearly as required. 

The proof of Theorem |7.2| is identical, only now since we do not have the 
hypothesis that the flow is outward initially we must take the moments with respect 
to p' — A(f> which is only equal to p outside of SI unless real analytic is assumed 
(as in Proposition |3.6[ l 

8 The Hele-Shaw flow and the complex HMAE 

In this section we will describe a connection between the Hele-Shaw flow and a 
certain complex Homogeneous Monge- Ampere Equation (HMAE). We will use 
this connection in order to prove short-time existence of the Hele-Shaw flow with 
varying permeability with empty initial condition. We recall that this extends a 
theorem of Hedenmalm-Shimorin |22| who prove this under the assumption k is 
real analytic and log k is strictly subharmonic. 

8.1 The Homogeneous Monge-Ampere Equation 

The complex Homogeneous Monge-Ampere Equation is a higher dimensional ana- 
logue of the Laplace equation on complex manifolds. We first describe some work 
of Donaldson (9) that gives an equivalence between regular solutions of a certain 
HMAE and families of holomorphic discs with boundaries in a so called "LS sub- 
manifold" which will be defined below. 

Let X be a compact Kahler manifold of complex dimension n, and let cjq be 
a Kahler form. This means that locally ujo is equal to idd4> where is a smooth 
strictly plurisubharmonic function. For example, if X has an ample line bundle L, 
and /i is a positive metric on L, then the curvature form iQh is a Kahler form on X 
which lies in the first Chem class of L. Any Kahler class in the same cohomology 
class as ujq can be written as cjq + idd(f) where is a global smooth function. 

Let D denote the closed unit disc and = D ~ {0}. Also let fio be the 
pullback of ujQ by the projection n: X x D — D. Suppose next that F is a 
smooth function on X x dD. We shall write Ft-(-) = F{-,t) thought of as a 
function on X, and suppose that uo + iddFr is a Kahler form for all r £ dD. A 
regular solution to the associated HMAE is an extension of F to a smooth function 
$ on X X D such that tJo + idd^r is a Kahler form for all r G D and 

on X X D. 

Now suppose $ is a regular solution to the above HMAE. Then kernel of 
r^o + idd^ defines a one dimensional complex distribution, which one can show 
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is integrable, thus giving a foliation of X x D by holomorpliic discs transverse to 
the X fibres (see [3J and |9|). 

A further connection with holomorphic discs comes about through the follow- 
ing idea due to Donaldson (9)- Let 8 := 0i + 162 denote a holomorphic 2-form 
on a complex manifold W, where Oi and 62 are real symplectic forms. A (real) 
submanifold V of VP' is called a LS siibmanifold if it is Lagrangian with respect 
to Gi while it is symplectic with respect to O2. 

We next sketch the construction of a holomorphic fibre bundle Wx over X 
such that Kahler forms in the same class as uJo are described as closed LS subman- 
ifolds in Wx- If 1^0 has a potential cf>u in some open set U, we identify Wu with 
the (1,0) part of the complex cotangent bundle. If Zi are local holomorphic co- 
ordinates any (1, 0)-form can be written as J]]^ ^idzi, thus (zi, ^j) are local holo- 
morphic coordinates of Wu- We also have non-degenerate holomorphic 2-form, 
namely O :— "^^i d^i ^ dzi. 

If V is some other open set where ixJq has a potential (fjv, then on U H V 
the transition function is given by 9(</f>v — 4'u)- It means that the section, locally 
defined as dffju, is in fact globally defined. By a simple calculation the graph of 
this section turns out to be a LS submanifold in Wx - Let oj be some other Kahler 
form in the same class as ujo- Then u — ujq + iddip for some smooth global 
function tjj, and the section d(j)u + dtp defines a section of Wx- By the same 
reasons this is also a LS submanifold. Donaldson shows that Kahler forms in the 
class of ixJQ are in one-to-one correspondence with the closed LS submanifolds of 
Wx- 

Now let _F be a smooth function on X x dD such that ujq + iddFg is positive 
for 6 G dD. For each 9 G dD, denote by Ae denote the corresponding LS 
submanifold in Wx - 

Theorem 8.1 (Donaldson (9)). There is a regular solution to the complex HMAE 
with boundary values given by F ijf there is a smooth family of holomorphic discs 
Qx '■ D ^ Wx parametrised by x £ X with 

1. Ttig^iQ)) = X e X- 

2. for each 6 G dD and each x £ X, 

gx{0) G As; 

3. for each r G D the map x i— >■ 7r((?2:(r)) is a dijfeomorphism from X to X- 

In fact for fixed t G D, the image of the map x i~> gx{T) is the LS submanifold 
At associated to the Kahler form uJo + idd^r- 

Using this theorem Donaldson applies the theory of the moduli spaces for em- 
bedded holomorphic discs with boundaries in a totally real submanifold to deduce 
the following: 

Theorem 8.2 (Donaldson O). The set of boundary conditions F for which the 
HMAE above has a regular solution ^ is open in the topology. 

8.2 Short-time existence of the Hele-Shaw flow with 
empty initial condition 

We will now consider a HMAE as above with X — f'^ and cjq ~ ^fs, the Fubini- 
Study form on P^. If we in the usual way think of as C U {00}, then on C the 
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Fubini-Study form has local potential 0fs = log(l + \z\^)- The local potential of 
a Kahler form in the same class as ujfs is a smooth strictly subharmonic (j> such 
that (f> — log(l + \ z\'^) has a smooth extension to the point at infinity. 

Now fix such a potential with A(j> — 1/k, where «: is a positive smooth 
function that encodes the permeability for the Hele-Shaw flow. For any 9 £ dD 
set 

F{z,e) ~(t,{ez)-ci>Fs{z). (13) 

Under the natural identification of Wc with fibred over C, the LS manifold 
associated to 4'{9z) is 

Ae = {{z, ^^{Oz)) : 2 G C} = {{z,e^{ez)) : 2 € C}. 

Now consider the C* -action on given by 

p{\) : {z,w) ^ {Xz,\'^w), A e C* 

which also restricts to an action of 5*^ = dD. Thus p{9) for 6 £ maps Ae to 
Ai. 

We will prove below that the HMAE with boundary condition F has a regular 
solution. Assuming this for now, we show how to produce the desired solution to 
the Hele-Shaw flow. By Theorem |8.1 [ there is an associated family of holomorphic 
discs gx in Wfi. Since the function F{z,d) is invariant under the S'^-action p, 
it implies that the solution <I> is invariant, and in particular the restriction $o is 
5^ -invariant. Set cpo := $0 + log(l + \z\^). Then (f>o is an S^-invariant strictly 
subharmonic function, which implies that ^?^(0) = and 7^(2) 7^ when- 
ever 2^0. This means that go{0) ~ (0,0), while for any a; 7^ 0, the second 
component of gx{0) is not 0. 

Furthermore the family of discs gx are be invariant under p, i.e. 

p(9)gx {Ot) = gex (r) foi 9 edD,r e D. 

In particular 

goie) = p{9-'){go)il) foi 6 & dD (14) 

and since n{go{0)) = we see ngo must be the constant map to 0. So, by conti- 
nuity of the family g, we have that as long as is sufficiently small, the image of 
gx:D^Wri,xj^O, lies in Wc = 

So fix such an x with x =^ 0. Then by twisting the disc using the C* -action p 
we get a new punctured disc fx:D^-^Cx¥^, 

Mr) := p{r){gx{r)). 

Since the second coordinate of gx{0) is not zero fx extends to a map D : C x¥^. 
Moreover for 6 £ dD and t £ D we have 

fe.ir) = pir)gox{r) = p{re)gxier) = fx{er) 

so the holomorphic disc defined by fx depends only on t — \x\, which we shall 
now denote by Et. Observe that fx{0) = (0, co) so St meets {0} x at (0, 00). 
Moreover if r G dD then gxif) G A^ so the boundary of Et lies in Ai. 
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Recall that the map 7(2;) := ■Kgx{'i) is a diffeomorphism of P'^ to itself. Since, 
by invariance nfx{0) = nggx{i-) = 'y(Sx) for 6 G dD we deduce 

dQt = {TT.f4e) ■.eedD} = {'y(ex) -.eedo} 

which is precisely the image of the circle of radius t = \x\ under 7. In particu- 
lar this shows that, for small t > Q,Q,t is an increasing family of smooth Jordan 
domains containing the origin, and tt: St — >■ vr(Et) is an isomorphism and the 
limit Q,t as t tends to zero is the origin. Moreover by Theorem |3.10| the family 
Q.t ~ 7i"(St) satisfies the Hele-Shaw flow as required, after possible reparameteri- 
zation of t. 

Thus it remains to prove that there is regular solution to the HMAE with 
boundary data F in (|T3j. Observe that if (p happens to be rotation invariant, i.e. 
(I>{tz) = (p{z) for all r £ dD, the HMAE has a trivial regular solution, namely 
$(2, r) := (j){z) — log(l + \z\^). (This reflects the trivial observation that if the 
permeability k is radially symmetric then the Hele-Shaw flow with empty initial 
condition will be the trivial solution in which Qt is a disc centered at the origin.) 

Let 1/) denote a rotation invariant potential that we also assume to be equal to 
|2p in a neighbourhood of the origin. By Donaldson's theorem we know that there 
is a C^-ball of radius e > centred at <^o — log(l + \z\^) consisting of boundary 
conditions for which the HMAE has a regular solution. 

Now an arbitrary smooth strictly subharmonic potential <j) can be written as 

4>{z) = Re(a + bz + cz + dz^ + ez'^) + f\z\^ + o(|z|^), 

where / is some real positive number. Our requirement on (fi was that = 
so without loss of generality we can assume that this first term is zero since it is 
harmonic, and also by rescaling that / = 1. This means that on some small disc 
centered at the origin (l> lies in the C^-ball of radius e/2 around ip- 

Pick Si and ^2 to be small positive numbers so that (f} is e/2 close to (1 + 
Si)tp — S2, and so that on the boundary of the small disc (j) < {1 + 5i)tp ~ 82. 
Still, in a neighbourhood of 0, we will have that <j) > {1 + Si)tp — 82- Let fg be 
a smooth non-negative function on the real line of integral one and with compact 
support in a (5-neighbourhood of 0. Then let 4>s denote the subharmonic function 




max((^, (1 + 5i)ip -52 + A)/(A)dA. 



This new potential will be smooth and strictly subharmonic. If 5 is chosen small 
enough, (f>s = (f)'ma. neighbourhood of the origin, while in a neighbourhood of the 
boundary of the disc it will be equal to (1 + + c where c is some constant 
(the integral of A/A minus 52)- We can thus extend it by letting it be equal to (1 + 
Ji)')/' + c outside of the disc. For sufficiently small 5, 04 — (1 + cSi)?/; — c will have 
-norm of less than e , and thus the same is true for F : = ^-^^^^^ — log ( 1 + 1 2 ^ ) . 

This implies that the HMAE with F as boundary condition has a regular so- 
lution By scaling, letting (1 + 5i)ujo be the new background Kahler form, the 
HMAE with boundary condition F = (1 + 5i)_F will have (1 + (Si)"l> as a regular 
solution as well, completing the proof. 
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